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Abstract 

The Eigenvalue Moment Method (EMM), Handy (2001), Handy and Wang 
(2001)) is apphed to the Ha = + iX^ + iaX Hamiltonian, enabling the 
algebraic/numerical generation of converging bounds to the complex energies 
of the states, as argued (through asymptotic methods) by Delabaere and 
Trinh (J. Phys. A: Math. Gen. 33 8771 (2000)). The robustness of the 
formalism, and its computational implementation, suggest that the present 
nonnegativity formulation implicitly contains the key algebraic relations by 
which to prove Bessis' conjecture that the eigenenergies of the Hq Hamiltonian 
are real. The required algebraic analysis of the EMM procedure pertaining 
to this problem will be presented in a forthcoming work. 
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I. INTRODUCTION 



A. General Overview 

There has been much speculation on the mechanism responsible for symmetry breaking 
within a special class of invariant Hamiltonians. It has been argued by Bender and 
Boettcher (1998), based on a conjecture by D. Bessis, that the class of Hamiltonians of 
the form + (zX)" admits bound states, within the complex plane, with real discrete 
spectra. Their arguments show that the VT invariance of the Hamiltonian is reflected in 
the wavefunction, '^*{—x) — ^(x), resulting in real spectra. However, Delabaere and Trinh 

(2000) have emphasized that VT invariance of the Hamiltonian is not sufficient to prevent 
symmetry breaking solutions. For instance, the Hamiltonian P'^+iX^+iaX admits bounded 
{LP') solutions on the real axis, which can have complex energies for a < acruicai < (thereby 
breaking VT invariance), or real energies, for a > CKcriUcai- 

An understanding of the underlying mechanism for symmetry breaking has remained 
elusive, despite the numerous investigations on the above, and related problems, by Bender, 
Boettcher, and Meisinger (1999), Bender et al (1999), Bender et al (2000), Bender and 
Wang (2001), Cahceti (2000), Delabaere and Pham (1998), Handy (2001), Handy and Wang 

(2001) , Levai and Znojil (2000), Mezincescu (2000,2001), Shin (2000), and Znojil (2000), in 
addition to those already cited. However, the recent work by Dorey, Dunning, and Tateo 
(2001) presents one possible explanation. 

Our objective is to seek alternative (and less analytical) arguments that can possibly 
shed some light on this matter. In this regard, we have attempted to implement a novel 
"positivity quantization" formalism based on the recent works by Handy (2001) and Handy 
and Wang (2001). These in turn make use of the Eigenvalue Moment Method (EMM) 
originally developed by Handy and Bessis (1985) and Handy et al (1988a,b). 

Our results, as communicated here, are very impressive. We are able to generate con- 
verging bounds for the (complex) discrete state energies, and arbitrary a. We are able to 
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confirm the general results derived by Delabaere and Trinh, although our methods do not 
rely on asymptotic estimates, as theirs do. 

It is important to emphasize that the EMM approach generates an infinite hierarchy of 
(closed) algebraic inequalities. We investigate the consequences of these relations from a 
numerical perspective. This is an important first step in identifying the algebraic relations 
responsible for symmetry breaking solutions, within our formalism. Our results strongly sug- 
gest that a careful algebraic analysis of the underlying EMM relations will serve to identify 
the theoretical structure leading to a proof of Bessis' conjecture. This will be presented in a 
future communication. We emphasize that unlike numerical integration schemes, which do 
not necessarily provide an undestanding of the underlying theoretical structure of a problem, 
the algebraic/numerical structure of EMM can. 

Beyond this, the ability to generate converging bounds to complex energy levels is a 
remarkable feat in its own right. This motivates the present communication. In a related 
work by Handy, Khan, Wang, and Tymczak (HKWT, 2001), they show how the Multiscale 
Reference Function formulation (Tymczak et al (1998a,b)) can easily generate the (com- 
plex) discrete state energies, for arbitrary a. Their estimation methods generate energies 
that fall within the bounds given here. As such, the present bounding theory provides a 
confidence test for the, numerically faster, MRF method. Because of this, we only gen- 
erate energy bounds for important a values near the first complex- real bifurcation point, 
Oicriticai — —2.6118094, as predicted by MRF. 
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B. Technical Overview 



The recent work by Handy (2001) introduced a new formalism that transforms the one 
dimensional Schrodinger equation into a fourth order, linear differential equation for S{x) = 
regardless of the (complex) nature of the potential, V{x) = Vr{x) + iVi{x): 

where S^'^\x) = d^S{x), and E = Eu + iEj, etc. Through this fourth order equation, one is 
able to transform the quantization problem into a nonnegativity representation suitable for 
a Moment Problem (Shohat and Tamarkin (1963)) based analysis, utilizing the Hamburger 
moments 

/+0O 
dx xPS{x). (2) 
-oo 

Such methods were originally developed by Handy and Bessis (1985), and Handy et 
al (1988a,b), and used to generate rapidly converging bounds to the bosonic ground state 
energy of singular perturbation/strong coupling (multidimensional) systems. This moment 
problem based, "positivity quantization" , approach is generically referred to as the Eigen- 
value Moment Method. An efficient algorithmic implementation requires the use of linear 
programming (Chvatal (1983)). 

We emphasize that the only bounded (i.e. functions within the ^-representation, 
or L^, within the S'-representation) and nonnegative solutions to Eq.(l) are the physical 
solutions (Handy (2001)). This is because if ^i;e{x) and ^2;e{x) are independent solutions 
of the Schrodinger equation, for arbitrary E (and thus unbounded, except for the physical 
solutions), then |vI'i.E(a;)p, |\E'2;e(x)P, ^'i.£(a;)^'2;E(x), and ^'i;e(x)^'2.b(x), are the inde- 
pendent solutions to Eq.(l), assuming V is complex. If V is real, then all the solutions 
are real, and only the first three configurations are independent; leading to a third order 
linear differential equation for S (Handy (1987a,b)). Thus, because of the uniqueness of 
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nonnegative and bounded solutions, application of EMM to the relevant moment equations 
(i.e. Eqs.(47 & 51) for Eq.(l), or the alternate moment formulation discussed in Sec. II) 
will generate converging bounds to the physical energies (Handy and Bessis (1985), Handy 
et al (1988a,b)). 

Despite this, we discover that for problems with complex energies, Ej ^ 0, the derivation 
of, as well as the actual, moment equation obtained from Eq.(l), are not the most efficient. 
A more efficient generation of the required tt-moment equation can be derived by working 
within a broader framework involving three coupled differential equations for the probabil- 
ity density, S{x), the kinetic energy density function, P(x) = |^''(a;)p, and the probability 
current density, J{x) — ^(^)^ ^^^2^ ^^^^ ■ The first two expressions are nonnegative con- 
figurations. In particular cases, the probability current density will also be nonnegative. 
This is discussed in Sec. II. However, this is not the immediate focus of the present work. 
Instead, it is to show how the realization of a nonnegative, linear, differential representation 
for S leads to a very effective quantization procedure, capable of confirming the existence 
of symmetry breaking solutions. 

Working within the coupled system of equations for {5, P, J}, one can generate many 
more moment constraints for the Ej ^ case. It is in this sense that we say that the moment 
equation derived directly from Eq.(l) is "incomplete", when Ej 7^ 0. Generating as many 
moment constraints as possible speeds up the convergence rate of the bounds. However, as 
noted above, even if we work with the reduced set of moment constraints generated from 
Eq.(l), the bounds generated will converge to the unique physical answer, just more slowly. 

Working directly with Eq.(l) leads to a complete set of moment equation constraints 
when Ei — 0. This is the case examined by Handy (2001), corresponding to the -|- (iX)^ 
problem (assuming real spectra. Handy (2001)). Application of EMM analysis in this case 
yielded impressive bounds for the first five energy levels (Handy (2001)), as well as complex 
rotated versions of the Hamiltonian (Handy and Wang (2001)). 

The sensitivity of the EMM procedure, as evidenced by the excellent nature of these 
bounds, strongly suggest that Bessis' conjecture is more likely a consequence of some un- 
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der lying algebraic identity, than a more subtle analytic constraint. 

Application of EMM to the more general case, Ej ^ 0, also yields very good bounds 
to the complex eigenenergies, as shown in this work. We have implemented our bounding 
analysis up to a relatively low moment order. Our objective has been to affirm the 5— EMM 
formalism's relevancy as an effective bounding theory. 
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II. DERIVATION OF THE 5-MOMENT EQUATION 



We derive Eq.(l) in a manner different from that presented by Handy (2001). As noted 
previously, the present formahsm is better suited for problems with Ej ^ 0. 
Denote the Schrodinger equation by 



where the normalized Hamiltonian is — —ed^ + V{x), involving a complex potential, V, 
and complex energy, E. We make explicit the kinetic energy "expansion" parameter, e. As 
recognized by Handy (1981), within a moments representation, kinetic energy expansions 
become analytic. This will be an important component of the present theory. The following 
discussion assumes that e is real and positive. 

It is readily apparent that if V*{—x) = V{x), and \l/(a;) is a bound state with complex 
energy, E, then ^*(— x) is another bound state solution with energy E*. Thus, complex 
roots come in complex conjugate pairs, for VT invariant Hamiltonians. 



A. The Coupled Differential Equations for \^{x)\\ |*'(a;)|2, and ^i^Wi'-^^^'i^Wi^) 

Define the quantities S{x) = \^{x)\^ P{x) = |*'(x)|2, and J{x) = _ 
The expressions Ei(a;) = '^*Hx'^{x) + c.c, Ai(x) = **iy^^'(x) - c.c, and E2(a:) = 
*'*i/a;^(x) + c.c, satisfy (Handy and Wang (2001)) 



i/,*(x) =£;*(x). 



(3) 



Ei(x) = €(2F(x) - S'\x)) + 2{Vr{x) - Er)S{x) = 0, 



(4) 



eP'{x) + {Vr{x) - Er)S'{x) - 2{Vi{x) - E:)J{x) = 0, 



(5) 



and 



|Ai(a;) = {Vi{x) - Ei)S{x) + ed,J{x) = 0. 



(6) 
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For the potential in question, V{x) — ix^ + iax, if S{x), P{x), J{x) form a solution set 
corresponding to energy {Eji,Ei), then S{—x),P{—x),J{—x) is a solution set for energy 
{Er, —Ej). 

In obtaining Eq.(l) (i.e. e = 1), we differentiate Eq.(4) and use Eq.(5) to substitute for 
P'. Upon dividing the resulting expression by — Ej, and differentiating, we use Eq.(6) to 
substitute for J'. For future reference, we make this process explicit. 

All of the physical configurations {S*, P, J} are implicitly assumed to be bounded and 
vanish at infinity. Thus, 

J{x) = -{ed,)-'(^{Vj{x) - Ej)S{x)), (7) 
where = J^^dxy. In addition, 

P[S; x] = ^S"{x) - -^{Vr{x) - Er)S{x). (8) 

We then obtain: 

- e''d,P[S{x);x] + e{VR{x) - Er)S'{x) + 2{Vi{x) - Ei){d,)-' [{Vi{x) - Ei)S{x)) = 0. (9) 

In order to transform this into Eq.(l), we simply apply the operator dx y^^J^_^^ , resulting in 
a fourth order, linear, homogeneous, diffential equation. 

B. The S- Moment Equation 

Let us denote by {up, Vp, Wp} the Hamburger moments of the three functions 5", P, J, 
respectively. Thus, Up = dx x^S^x), etc., for p > 0. For the present problem, Vr = 0, 
and Vj — x^ + ax. Multiplying each of the three equations by x^, and integrating by parts, 
yields 

2evp - p{p - l)eUp-2 - 2ErUp = 0, (10) 
pevp^i + ErpUp^i - 2{wp+s + aWp+i - EjWp) = 0, (11) 
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and 

+ aup+i - EiUp) - pewp^i = 0, (12) 

p > 0. Note that from Eq.(lO) (for p = 0) must be positive. 

We can convert this into a moment equation for the {up}'s by using the first moment 
relation to solve for Vp in terms of the it's. Likewise, taking p — > p + 1 in the last moment 
relation, determines Wp in terms of the u's. Finally, substituting both relations in the second 
equation generates a moment equation for the u^s: 

2/ p3 ^ ^ ^ 2Ej ^ / 1 1 \ 

e (P - ^ + -j)Up-3 + 2peERu{p - 1) - j^^Up+i + 2aEi[-^ + -—^jup+2 - 

p + ^p + 2 p + A'' p + A 

This moment equation holds for all e, including e = 0. In the latter case, upon multiplying 
the M- moment equation by {p + l){p + 2){p + A), the resulting relation (quadratic in p) 
incorporates the relation given in Eq.(12) for e = 0. 

We note that the it-moment equation does not include one important additional moment 
constraint, that for the p — relation in Eq.(12). This yields 

^3 + aui — EjUq. (14) 

Of course, this is the relation one obtains directly from the Schrodinger equation, upon 
multiplying it by ^*{x): 

J dx (eP(x) + V(x)S(x)) ^Ejdx S(x), (15) 

and identifying the real and imaginary parts (the latter corresponding to Eq.(14)). 

The recursive, linear, homogeneous, structure of the ti-moment equation tells us that all 
of the moments are linearly dependent on the first seven moments {m^IO < i < 6}. The 
additional constraint in Eq.(14) allows us to solve for in terms of Uq and Ui. Thus, the 
reduced set of independent moments is {ue\0 < £ < 2,4 < £ < 6}. These are referred to as 
the missing m,om,ents. In addition, we must impose a suitable normalization condition. The 
details are given in the next section. 
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There is an important theoretical point that must be stressed. The recursion relation in 
Eq.(13) is the moment equation resulting from integrating both sides of Eq.(9) by x^. Both 
Eq.(9) and Eq.(l) are equivalent to each other, as explained above. However, the moment 
equation generated by Eq.(l) (refer to Eq.(47) and Eq.(51)) is different from Eq.(13), as 
explained in the Appendix. It will involve one more degree of freedom (i.e. missing moment 
order) than Eq.(13). However, since Eq.(l) and Eq.(9) are equivalent, and Handy (2001) 
has argued that the only bounded and nonnegative solutions are those corresponding to the 
physical configurations, apphcation of EMM to either will result in converging bounds. Nat- 
urally, the moment equation involving fewer independent variables (i.e. missing moments) 
will yield faster converging bounds. 

In addition, because Eq.(13) is the moment equation of Eq.(9), and similarly for Eqs.(47 
& 51) and Eq.(l), any additional constraints, such as Eq.(14), are not required in order 
to generate converging bounds. Such constraints only improve the convergence rate of the 
bounds by reducing the number of independent, missing moment, variables. 

^From a different perspective, the manifest difference between Eq.(13) and Eq.(47 & 51) 
hinges on the fact that in obtaining Eq.(l) we implicitly take p — > 1 in Eq.(ll). We note 
that this implies that Eq.(l) cannot generate the extra constraint + awi — EjWq — 0, 
that is, the p = relation from Eq.(13). 

If Ej — 0, and the nonnegative configurations S{x) and P{x) are symmetric, then these 
extra constraints (including Eq.(14)) are nonexistent. That is, Eq.(l) yields a complete set 
of moment constraints if £'/ = 0. 

The moment equations in Eqs. (13-14) are the preferred relations ( because they are easier 
to derive, and involve less missing moments), if Ei ^ 0. 

We will work with the w-moment equation as given above, complemented by Eq.(14). 



10 



C. The Zeroth Order e-Contribution 



As previously noted, one of the most important reasons for working within a moments' 
representation is that it is analytic in e. 

li e ^ 0, there is equivalency between the sets of equations Eq. (10-12) and Eq. (13-14). 
When e = (and thus Eji — 0), the w moments decouple from the u moments. We find 
that Eq.(13-14) includes more solutions than those generated from Eq.(12). However, not 
all of these will be consistent with EMM quantization. One formal way of understanding 
this is to consider Eq.(9) when e = 0. We obtain 

{Vj{x) - Ei){d,)-'{{Vj{x) - Ej)S^'\x)) = 0. (16) 

For real, bounded, configurations, {f{x),g{x)}, the integral dx f {x){dx)~^ g{x) 
becomes dx (((5a;)"^)"''/(a:)^ di^), where {{dx)~^y = f^°° dx^. This follows from 
/_+~ dx fix) dx, g{x,) = /_+^ dx d^ fix) g{^ + x) = /_+^ dx di g{x)f{x-0 = 
f^^dx g{x)!:^dx, fix,). 

Multiply Eq.(16) by x^Rfs^x), where Rpix) is a regulating (bounded) function which 
reduces to unity when /3 ^ 0. Integrating, and using {(^dx)~^^\ gives (i.e. in the ^ 
limit) 

((9J-^)tfa;f+^ + axP+^ - Ejxp) ^ -^xP+^ + -^xP+^ - -^x^+K (17) 

This, in turn, upon multiplying by {Vi{x) — Ei)S^^\x), and completing the x- integration, 
gives the zeroth order (in e) moment relation in Eq.(13). 

Therefore, the e = moment equation in Eq.(13) (i.e. corresponding to Eq.(16)) tells us 
that 

{Vi{x)-Ej)J^'\x)^0. (18) 

This has the general, formal, solution J^^\x) — J2e JiS{x — Ti), where the r^'s are the turning 
points: 
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Vi{n{E)) = Ei. (19) 

The J'/s are arbitrary. However, the e = solution to Eq.(12) really correspond to J = 0. 
We can (formally) argue this by applying dx{Vi{x) — Ej)"^ to Eq.(18), yielding 

d,J^'\x) = {Vi{x) - Er)S^'\x) = 0, (20) 

which is the underlying configuration space relation corresponding to Eq.(12), for e = 0. 
The only possible, bounded, solution is J = 0. 



Thus, the only solution to Eq.(13), for e = 0, consistent with the EMM quantization 
constraints (which demand boundedness and nonnegativity) , as discussed in the following 
section, should be that corresponding to J^^^ (x) — 0. 

Consistent with the previous discussion, the e = moment equations (Eq.(10-12)) yield 

+ - Eiw^'^ = 0' (21) 

and 

p > 0. The solution set to these are Up — ^(,Ae,Tf, and Wp — YliiJi'^^- 

All of the (complex) turning points, for arbitrary Ej^ contribute to the zeroth order 
structure of S{x). 

D. The Moment Problem Constraints 

The Moment Problem conditions for nonnegativity (Shohat and Tamarkin (1963)) cor- 
respond to the inequalities: 

N 

dx {^CnX^fS{x)>0, (23) 

foTN<oo and C„ arbitrary, real, variables. These are usually transformed into nonlinear 
(in the moments) Hankel-Hadamard (HH) determinantal inequalities, Aq^n{u) > 0, where 

12 



^o,Ar(w) = Det{ui+j), and < i, j < N. However, we prefer their linear (in the moments) 
equivalent, corresponding to the quadratic form inequalities: 

N 

Cn, > 0. (24) 

ni,n2=0 

These will be referred to as the linear HH nonnegativity constraints. 

The moment ("Hankel") matrix Um+n^ is linearly dependent on the missing moments, 
{uo,ui,U2,U4,U5,Uq}, nouliuearly dependent on {Eii,Ej), and analytic with respect to e. 
With respect to the following discussion, we note that through an appropriate normaliza- 
tion prescription (as discussed in Sec. Ill), the missing moment variables will lie within a 
bounded, convex, domain. 

The physical solution must satisfy all of the above constraints. To any finite order. A'', 
if at a given energy value, E — {Er, Ej), there is a missing moment solution set, Un-e, then 
it must be convex. The EMM eigenenergy bounding procedure simply involves determining 
the energy subregions, {Er, Ej) e TZjq-j, for which Un;e exists. The j-index enumerates the 
discrete states. 

Based on the many apphcations of EMM over the last sixteen years, one expects the TlN-j 
regions to be connected and bounded, although not necessarily convex. This is supported 
by our empirical results. 

The boundary of the smallest rectangle containing TZ^-j (i.e. [E^^\e^I^''] x [E')["\Ef^] D 
TZn-j) define the bounds for {Er, Ej) e TZn-j- In practice, we numerically determine a slightly 
larger rectangle than this. Our numerical analysis will be based on the above relations. The 
following discussion addresses an interesting side issue that enhances our understanding of 
the above relations. 

The numerically minded reader may wish to skip to Sec. III. 

E. Simplification of the EMM Constraints for Nonnegative J{x) 

The Hankel-Hadamard inequalities are automatically satisfied by the atomic distribution 
S{x) 5(°)(x) = j:iAeS{x - n), provided At > 0. The zero equality is satisfied by 
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polynomials, Pn;c{x) = Z^^=o CnX" whose roots include all of the turning points, Pn;c{x) — 



For physical (non-atomic distribution) solutions, only the strict inequality can be satis- 
fied. In this case, we can work with the relations 



Normally, one would prefer to work with such positive relations because they will not contain 
any zeroth order e dependence; thereby generating the positivity constraints that really con- 
tribute to quantization. However, in the present case, such inequalities are not independent 
of zeroth order e contributions. 

Instead, as argued before, the zeroth order structure of the moment equation in Eq.(13) is 
due to the probability current, J. However, we can only work with nonnegativity constraints 
for the current, if it is nonnegative. 

If the probability current is nonnegative, then 



has no zeroth order e dependence. That is, from Eq.(18), all the zeroth order terms are 
elliminated. This set of constraints is algebraically simpler than working with the analogue 
of Eq.(23), as apphed to J, assuming J > 0. 

F. Properties of the S, P, J Equations and (Minimal) Conditions for J's Positivity 

We assume that e > 0, Eji > 0, and Vr — 0. The turning points satisfy V/(t^) — Ej. The 
following analysis is restricted to the real axis. We also assume that all of the configurations 
have analytic extensions into the complex-x axis. Thus, S{x) = \E'(a;)\l'*(x*) is the analytic 
extension of |^(a;)p, etc. Since V{x) is a regular function, the wavefunctions ^(x) and 
^*{x*) are regular functions in the complex x-plane; hence, so too are s{x), P{x), and J{x). 

Both S and P are nonnegative. Because of their definitions, both S and P cannot be zero 
simultaneously, except at infinity (i.e. otherwise ^ and ^' would be zero simultaneously. 



PM-,-d^)^iyi{x)-Ej). 




(25) 




i=o 



(26) 
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generating the trivial zero solution). When S{xo) — 0, or P{xo) — 0, then "^{xo) — 0, or 
^'(a^o) = 0, respectively; therefore, the probability current, J, is zero at all of the zeroes of 
S and P. 

It follows from the nonnegativity of S that if S{xo) — 0, then these are extremal points, 
S'{xo) = 0. The same holds for P. 

Lemma # 1: If S{xo) = 0, and Xo 7^ r ( a real turning point), then J becomes 
negative in the neighborhood of Xo 

Since S{xo) = 0, then S'{xo) = 0, and J{xo) = 0. ^Prom Eq.(6), J'{xo) = 0. Dif- 
ferentiating Eq.(6) once, we obtain J"{xo) = 0. Differentiating it a second time yields 
-eJ"'{xo) = {Vi{xo) - Ei)S"{xo). ^From Eq.(4) 2P{xo) = S"{xo) ^ 0. Thus, since 
Xo 7^ r, it follows that J"'{xo) 7^ 0. That is, the local power series expansion for J be- 
comes J{x) = \j"'{Xo){x — XoY + 0{{x — XoY). O 

Lemma # 2: If J'{xo) — 0, and Xg 7^ r, then J bececomes negative in the neigh- 
borhood of Xo 

^Prom Eq.(6) it follows that S{xo) — and Lemma #1 applies, o 

Lemma # 3: If {Vi{x) — Ej) is Eisymptotically monotonically increasing, with one 
real turning point, then J{x) > 0. J{x) is strictly increasing for x < r, and strictly 
decreasing for x > r. 

By assumption, there is only one real turning point, r. Also, \iva.x^±^{Vi[x) — Ej) — 
±00. ^Prom Eq.(5) we have (eP + ErS)' — —2{Vi{x) — Ei)J{x)] however, since S and 
P are asymptotically positive and decreasing to zero, we have lima;^±oo(c-P + E^S)' = 0"^. 
Accordingly, limx-^±^ J{x) =0"*". It then follows that there must be two points (possibly 
the same), coming in from —00 and -|-oo, where J is a positive local maximum. Denote 
these by Xg^ < Xo2- At these points we have J'{xoi2) — 0- If either of these points is not r, 
then by Lemma #2, we have a contradiction. Thus 2 — 'T-i J{x) is strictly positive. 
Clearly, J'{x) > 0, for x < r, and J'{x) < 0, for x > r. o 



15 



Lemma # 4: If + f £;| > 0, then J{x) > 

Assume V/(.r) = -\- ax. If a > 0, then Vi{x) — Ej satisfies the conditions of Lemma 
#3. If q; < 0, then Vi{x) will have a local maximum at a:_ = — (^)^, and a local minimum 
at x^ — —X-. At these locations, we have Vi{x-) — and V/(x+) = — yf(x_). If 

Ei > Vi{x-), or Ei < Vi{x+), then the conditions of Lemma # 3 are satisfied, o 

Lemma # 5: If J{x) > 0, then its local extrema must occur at the turning points 

The local extrema correspond to J'{xo) = 0. U Xg ^ t (where there can be more than 
one real turning point), then according to Lemma # 2 we contradict the assumption that 
J{x) > 0. Thus, for nonnegative J's, its zeroes must coincide with the turning points, o 

Lemma #6: eP < E^S within the local mcixima regions of S, and eP > E^S 
within the local minima regions of S. They intersect at 5''s inflection points. 

This immediately follows from Eq.(4), eP = ErS+^S". We see that eP{x) and ErS{x) 
intersect at the infiection points, S"{xi) — 0. At infinity, since Yim.^^^^ S" {x) — 0+, then 
lim3.^-i-oo(€P(x) — ErS{x)) = 0"*". Between any two successive extremas of EjiS{x) there 
must be an intersection by eP{x). o 

Lemma # 7: Let (V} — Ei) be asymptotically monotonically increasing, with one 
real turning point (r). Let {xi,X2) be an interval whose endpoints correspond to 

successive extrema for S*, where S'{xi) = S'{x2) = 0. If {xi,X2) C (r, oo), and S{x) 
is increasing within the interval, then P{x) must be monotonically decreasing 
within the interval, and vice versa (i.e. S <-> P). When {xi,X2) C (— oo,t), then if 
S is decreasing within the interval, P must be monotonically increasing within 
the same interval, and vice versa. 

Under the conditions of the Lemma, J is positive (by Lemma # 3); therefore, at any 
extremal value for S (i.e. S'{xo) = 0) , from Eq.(5), it follows that eP\xo) = — 2(V/(.To) — 
Ei)J{xo). If Xo > T, then {Vi{xo) — Ej) > 0, and we have that P'{xo) < 0. Likewise, for S'. 
That is if Xo > r, and P'{xo) = 0, then from Eq.(5) ErS\xo) = -2{Vi{xo) - Ei)J{xo) < 0. 
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Let T < xi < X2, and let {xi,X2) denote an open interval defined by two successive 
extremum points for S, within which S is increasing. It then follows that P'{x) must 
be negative at each endpoint. However, P cannot have a local extremum within the 
interval,P'(a;e) = 0, since then S'{xe) < 0, contradicting the assumptions made. Thus, 
P must be stricly decreasing within the closed interval [xi , X2] ■ 

If xi < X2 < T, then if {xi,X2) denotes an open interval on which S is decreasing 
{xi^2 being two successive extremum points), then P'{x) must be positive at each endpoint. 
P' must remain positive within the interval, otherwise at any internal extremum point, 
S'{xe) > 0, contradicting the assumptions made, o 

^Prom Eq.(5), under the conditions of Lemma # 3, then eP{x) + EjiS{x) has only one 
extremum point, a global maximum, at the turning point. 

Lemma # 8: If P{xo) — 0, then S{x) has a local mciximum at Xo 

If P{Xo) = 0, then P'{xo) = and J{Xo) = 0. From Eq.(5), S'{Xo) = 0. From Eq.(4) we 
have that eS"{xo) + 2EjiS{xo) — 0. This cannot be satisfied at any local minimum (since 
there, one has S" > 0, and S >0). We cannot have S" — and -S" = 0, since both S and P 
cannot be simultaneously 0. The only possibility is that Xg is a local maximum for S o. 

Lemma 7^ 9: If S{xo) = 0, then P{x) has a local maximum at Xo 

If S{xo) = 0, then S'{xo) = and J{xo) = 0. From Eq.(5), P'{xo) = 0. ^From Eq.(6) 

J'{xo) = 0. If we differentiate Eq.(5), eP"{xo) = -ErS"{xo). From Eq.(4), 2P{xo) = 

S"{xo) > 0, therefore P"{xo) < 0, corresponding to P having a local maximum at Xg. o. 
The preceding Lemmas do not shed any immediate resolution to the question as to when 

is Ej — 0. This is because they are mostly of a local nature and do not make use of the 

boundedness criteria, for physical solutions, other than Lemma's #3 and #4. Nevertheless, 

these represent two important contributions. 

A potentially significant question is, for what {a, En, Ei) values will the quantization of 

J > 0, through EMM, yield any results? The impact of J's nonnegativity on the reality of 

E is presently under investigation. 
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For completeness, we note that if Ei — 0, and Vi is asymptotically monotonic, with one 
zero point, then J > 0, as is evident from Eq.(6). 



G. Additional Properties 



Let V{x) be an arbitrary, VT invariant, potential: V*{—x) — V{x). For any bounded, 
complex, solution to the Schrodinger equation, \l/(a;), with energy, the configuration 
^*(— x) solves the same quantum problem but with eigenenergy E* . This means that 
both {S{x),P{x),J{x)} and {S{-x),P{-x), J{-x)} solve Eqs.(10-12) for {Er,Ei) and 
{Er, —Ei), respectively. That is {up,Vp,Wp} {(— l)^'Mp, (— (— satisfy the 
coupled moment equations (together with {Eji,Ei) {Er, —Ei)). 

There is another perspective on the above. The ^*{x) solves the equation 



Thus, we may regard as the solution to the Schrodinger equation upon taking e — > — e 

(with eigenenergy —E*). Similarly, S{x), P{x), —J{x) solves Eqs. (10-12) for e — e, Er 
-Er, and Ei Ei. 

The Eigenvalue Moment Method is an quantization theory which is analytic in e. That 
is, bounded configurations are normalized according to integral relations (i.e. / dx S{x) — 
1, etc.) and not the usual conditions within the usual quantum mechanics Hilbert space 
formulation. In the present case, since we are working with {S, P, J}, our normalization 
coincides with the formulation of quantum mechanics, so long as e 7^ 0. Within the 
usual quantum mechanical formalism, the eigenstates of the position operator (i.e. the 
translated Dirac delta function) are non-normalizable (within the norm). However, they 
are normalizable within the norm inherent to the EMM approach. 

The final observation is that upon performing the change of variables ?/ = 7, the Hamil- 
tonian under consideration becomes 



edl<i/*{x) + V{x)'i/*{x) 



E*^*{x). 



(27) 



-a^* + i{s^y^ + asy)"^ = E^ . 



(28) 



18 



Dividing by s^, we see that 



i;(4,4) = ^. (29) 



So it follows that if s = es , then 



3 d 

E{a,e)^e-^E{-,1). (30) 

65 
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III. EMM-NUMERICAL ANALYSIS OF THE + jx^ + alX HAMILTONIAN 



The recursive, linear structure, of the ti-moment equation in Eq.(13) can be written as 
(i.e. e = 1) 



(31) 



for p> 0, where Mp^i satisfies Eq.(13) with respect to the p index, as weU as the initiahzation 
conditions Me^^e^ = Se^^e2, for < £1,2 < 6. 

These relations are supplemented by Eq.(14), written in the form 1^3 = EjUq — aui. We 
can then substitute in the previous relation, obtaining 

6 



(32) 



where 



' Mp^o + EiMp^s, £ = 
Mp,i-aMp,3, e^l 



M„ 



P,2) 



e = 2 



(33) 



0, £ = 3 
[Mp,e, 4<£<6 

Since Mp^a = 0, we will work with the reduced set of independent moment variables 
uo — uq, v\ — vi — U2, i^z — M4, ^4, — U5, and — uq. Thus 

5 



(34) 



where flp^g — Mp^^, for < £ < 2, and flp^^ — Mp_^+i, for 3 < £ < 5. 

Finally, we define our normalization condition with respect to the even order Up moments 
J2i=o = 1- Since each of these is positive, it insures that the requisite linear programming 
analysis is done within the five dimensional unit-hypercube [0, 1]^. 

In order to impose this normalization condition, we invert the <-> U2(. relation (i.e. 
U2t = Et=o ^2i,e,i^i^, < £ < 5) 
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ev=o 



(35) 



< ^ < 5, and substitute into Eq.(34), Up = ELo ^p/(Et=o -^^A'^^s^J obtaining 



5 




where Tp^^ — ^4=0 We now insert the normahzation condition (i.e. solve for uq 



in terms of the first five even order moments) , obtaining 




for p > 0. The hnear programming EMM algorithm is implemented on these relations, 
within the context of Eq.(24), or the equivalent quadratic form counterpart to Eq.(25). 

The numerical results of our analysis are given in Table I, for a selected number of a val- 
ues, of interest within the asymptotic analysis by Delabaere and Trinh (2000). The reasons 
we do not quote more bounds, for more a values, is that the Multiscale Reference Function 
(MRF) analysis of Handy, Khan, Wang, and Tymczak (HKWT, 2001) is numerically faster, 
and yields results lying within the EMM bounds given here, for the selected a values. This 
strongly suggests that the MRF analysis is correct. This is a good example of how the 
present "bounding" method can be used to test other (generally much faster) estimation 
methods. 

The MRF approach predicts that, for the lowest lying discrete states, there is a critical 
a value below which complex energies appear. This is given by 



The data given in the Tables (particularly Table II) is meant to test the reliability of 
this. It is clear that it does. 

The bounds for the real and imaginary parts of the eigenenergies are very good. In the 
Tables, Pmax defines the maximum moment order used. That is, it is the total number of 
Hamburger moments used (i.e. {//p|0 < p < Pmax})- H the EMM procedure is applied to 



^critical = -2.6118094. 



(38) 



21 



the Ei — case (i.e. Handy (2001)), corresponding to symmetric 5'(a:)'s, then the formahsm 
converts to a Stieltjes moment representation which only involves the even order moments: 
{a*2p|0 < P < Pmax}- Thus, a Hamburger moment order of Pmax-, corresponds to a Stieltjes 
moment order of P^}^ — The tight bounds in Handy 's original work (2001) required 

p^^i-om. 
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TABLES 



TABLE I. Bounds for the Discrete States of + ix^ + iaX 



a 




Ef> < En < 4"' 


) < < 


-3 


20 


0.7 <Er< 1.7 


0.4 < ±Ei < 1.0 


-3 


24 


1.10 <Er< 1.45 


0.5 < ±Ei < 0.9 


-3 


28 


1.20 <Er< 1.23 


0.72 < ±Ei < 0.77 


-3 


32 


1.219 <Er< 1.230 


0.756 < ±Ei < 0.768 


-3 


36 


1.224 <Er< 1.228 


0.758 < ±Ei < 0.762 


-3 


40 


1.22561 <Er< 1.22608 


0.75980 < ±Ei < 0.76055 


-3 




1.225844* 


.760030* 


-2 


20 


.416 < < .719 


-.5 <Ei<.5 


-2 


24 


.607 <Er< .636 


-.03 < < .03 


-2 


28 


.610 <Er< .625 


-.5 X 10-2 <Ei <.5x 10- 


-2 


32 


.619 < < .625 


-.2 X 10-2 < £;/ < .2 X 10- 


-2 


36 


.6203 <Er< .6213 


-.45 X 10-=^ <Ei < .45 X 10- 


-2 


40 


.62083 < ^i? < .62105 


-10-^ <Ei < 10-"^ 


-2 




0.6209137* 


0* 



*Multiscale Reference Function formulation by Handy, Khan, Wang, and Tymczak (2001) 



23 



TABLE II. Bounds for the Discrete States of + iX^ + iaX 



a 


p 

^ max 




) <Ei< < ) 


-2.610 


20 


517 < En < I 920 


-.6 < 


< 


.6 


-2 filf) 


24 


Q4n <" Kd ^ 1 800 


-.4 < £;/ 


< 


.4 


-2.610 


28 


1 083 < E'd < 1 586 


-.2 < Ei 


< 


.2 


-2 610 




1 21 1 <r Ed < 1 361 


-.08 < 


< 


.08 


-2.610 


36 


1 214 < i?t? < 1 250 


-.23 X 10"^ < Ei 


< 


.23 X 10-1 


-2.610 


40 


1 2135 < En < 1 2617 1 2617 < ii^o < 1 3581 


-.5 X 10-2 < Ei 


< 


.5 X 10-2 


-2.610 


42 


1 2317 < £^7? < 1 2367 1 3179 < < 1 3356 


-.25 X 10-2 < Ei 


< 


.25 X 10-2 


-2.610 


42 


1 234216* and 1 332059* 


0* 






-2.614 


20 


515 < ii^D < 1 925 


-.525 < Ei 


< 


.525 


-2.614 


24 


953 < En < ^ 808 


-.39 < Ei 


< 


.39 


-2.614 


28 


1.083 < < 1.589 


-.12 < £;/ 


< 


.12 


-2.614 


32 


1.238 <Er< 1.326 


.01 < ±Ei 


< 


.11 


-2.614 


36 


1.256 < £;r < 1.309 


.030 < zbE"; 


< 


.065 


-2.614 


40 


1.278 <Er< 1.286 


.050 < ±Ei 


< 


.065 


-2.614 


40 


1.282333* 


.0538739* 





*Multiscale Reference Function formulation by Handy, Khan, Wang, and Tymczak (2001) 
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Our results (combined with those of HKWT (2001)) confirm the asymptotic estimates 
provided by Delabaere and Trinh. In particular, we can access regions in the a parameter 
space which were difficult within their formulation. 

The numerical implementation given, of the present formalism, is meant to suggest the 
power of the method. Our intention is not to present an exhaustive numerical analysis over 
a wide range of parameter values. 

A. Basic Algebraic Structure of the u-Moments 

Upon using Eq.(14) to solve for in terms of Mo,i, and incorporating this into the 
moment recursion equation in Eq.(13), we obtain 

/25a£;? ^ X 2ba^Ei 5Ej Sa^ TEj 8a 
Us = [ — g h 5eEiijuo — ui —U2 —U4 + - -^"e, (39) 



and 



uw=( ^ + ^ ~ UaeEnjuo + ( ^ - llE^jui 

/45aE'j \ / Q o\ 27aET o . v 

+ (— + 21eER)u2 + [ia^ + 12E])u4 r— + Sa^e- (40) 



(41) 



In obtaining Eq.(36), we must invert the relationship {us,uiq} {uijU^}. This takes 
on the form 

Us-R8[uo,U2,U4,U(,] \ / -25^^,Z|I \ / „^ 

U,o-Rlo[uo,U2,U4,Ue]) V _ 11^3 _ 27^ J [^^ 

where -^8,10 denote the remainder terms in the previous relations. The determinant of the 
above matrix is 

313^ ' WrE ll=2«^^l + ?^l- (42) 



Accordingly, 



U5, 



(43) 
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Thus, when Det{Ei) — 0, the even order moments {t^o, ^2, U4, uq, us, uiq} are constrained to 
satisfy two additional relations. That is, the effective dimension of the system drops from 
6 — > 5 (before imposing a normalization constraint). The impact of this on the EMM bounds 
is unclear, at the present time. However, a convenient feature is that Det{Ei) > 0, under 
the conditions of Lemma #4. 
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IV. CONCLUSION 



We have presented a nonnegativity representation formalism, amenable to the Eigenvalue 
Moment Method, for generating converging bounds to the (complex) discrete eigenenergies 
of one dimensional, VT- invariant Hamiltonians. Our analyis was presented within the 
specific context of the Ha. — P"^ -\- ix^ + ioLX Hamiltonian, previously analyzed by Delabaere 
and Trinh. Our formalism readily confirms (numerically) the existence of both symmetry 
breaking solutions, and symmetry invariant solutions. The preliminary results given have 
focused on the properties of the low lying states near the first symmetry breaking bifurcation 
point, with respect to the a parameter, as predicted by the MRF eigenenergy estimation 
method of Tymczak et al (1998a,b). Our preliminary, yet highly accurate, numerical results 
are consistent with the asymptotic analysis methods of Delabaere and Trinh. 
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VI. APPENDIX: EMM ANALYSIS OF EQ.(l) 



In this section we focus on the structure of Eq.(l) and the moment equation resulting from 
it. This is not the most efficient approach, as indicated earher. Instead, it is best to derive 
the moment equation for S by first working with the coupled S, P, J equations, deriving the 
corresponding coupled moment equations, and then reducing these to one moment equation 
for S. Nevertheless, for the sake of completeness, we outline the issues that need to be 
addressed if one wants to go from Eq.(l), directly, to a moment equation. 

The fourth order equation for the + ix^ + aix Hamiltonian is 

(44) 

where A(a;) — Ej — x{a + x^). 

The function coefficient introduces singularities on the physical domain, A(r) = 0. 
However, all of the solutions to Eq.(44) will be regular (Handy (2001)). These roots are the 
effective "turning points" of the differential equation. They satisfy 

+ ar = Ei. (45) 

The extremal points of the function Vi{x) = x^ -\- ax satisfy x^ — — |. 

If q; > 0, then there is only one real root to Eq.(45). If a < 0, then there could be one or 
three roots, depending on Ej. That is, x^^^ = ±y and Vi{x^^^) = =F/5, where {3 = 
If l-E'/l < or l^/l < .3849|q;|2, then there are three roots. 

It is best to translate Eq.(44) to any one of the real r roots, in order to simplify the 
structure of the ensuing moment equation. Therefore, we will work with the translated 
variable, ^ — x — t. Let us define the power moments 

I^P=[d^eS{0, (46) 

where C is a contour in the complex-(^ plane that sits on top of the real-^ axis and deviates 
around the origin. Since all the solutions (physical or not) to Eq.(44) are regular in ^, we see 
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that if 2? > 0, then for the physical solutions, we can deform the C contour to be identical to 
the real-^ axis. If p < 0, we cannot do this, and must retain C. The p > power moments 
are referred to as the Hamburger moments. 

We note that A(0 = C^(0> where T(0 = Sr^ - a + 3r^ + If there are three, real, 
r-roots, then T({) has two zeroes along the real axis. If there is only one r root, then T has 
no zeroes along the real axis. 

Multiplying Eq.(44) (translated by an amount r) by ^^(A(^))^, and integrating by parts 
over the contour C, we obtain a moment equation valid for — oo < p < +oo: 

- (3t' + a)p{A -Ap-p''+ p3)/Xp_3 - 3tp(-4 -p + Ap^+ p^)np-2 
-pip + 2)(6 + 12t^Er + AaER + 7p + p')^lp-l - 12tEr{A + 5p + p')iip - AEr{12 + 8p + p')iip+x 

+0//p+2 + 4(3t^ + afjip^^ + 36t(3t^ + afjip+i 
+12(36t^ + ISar^ + a^)/Up+5 + 36(9r^ + 2Ta)np+6 + (144t^ + 12a)np+7 + SGr^p+s + 4//p+9 = 0. 

As argued by Handy and Wang (2001), as well as Handy, Trallero, and Rodriguez (2001), 
in order for EMM to yield converging bounds, it is important that the proper (Hamburger) 
moment equation uniquely correspond to the desired system. Let us represent Eq.(44) as 
O^S{^) — 0. Now consider a more general problem corresponding to O^S{^) — T>(^), where 
the inhomogeneous term corresponds to a distribution like expression, supported at the 
zeroes of A(^). When we multiply this system by A(^)^, the corresponding inhomogenous 
term can (effectively) disappear. If one is not careful, and solely restricts the moment index 
to nonnegative values, then the resulting moment equation cannot distinguish between the 
desired problem (corresponding to "D = 0) and those for which V ^ 0. In such cases, the 
EMM algorithm will not generate any bounds. 

In order to insure that our moment equation refers to T? = 0, we must work with 
the moment equation evaluated for p > —2, which represents the most singular function 
coefficient remaining after Eq.(44) is multiplied by T(^). In the work by Handy and Wang 
(2001), the moment equation for p = —2, —1 yielded additional contraints for the Hamburger 
moments. The same is true in the present case. Thus, the moment equation for p — —2 
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becomes 



4(3t2 + afni + 36r(3r2 + affX2 + 12(36r^ + ISar^ + a^)n3 
+36(9r3 + 2tq!)//4 + {Uir^ + 12a)fX5 + 36t/X6 + 4/// = A*-4), 



(48) 



where 



E(//_2,Ai-4) 



(24T£;ijyLt_2 + 36r//_4). 



(49) 



The moment equation for p = —1 becomes 




(50) 



By combining the previous two relations, we can express /Us in terms of the Hamburger 
moment {/Uo, A*!, ■ ■ ■ , fJ-v}- 



This is an important additional constraint on the Hamburger moments, as noted in 
the work by Handy and Wang (2001). Since the only bounded and positive solutions to 
Eq.(44) are the physical solutions, EMM will work directly on the moment equation in 
Eq.(47) supplemented by the above relation (which is essential, otherwise no bounds will be 

generated). 

The EMM numerical implementation of the above would involve (for fixed a) using r 
as the variable parameter (and then computing Ej from Eq.(45)), in addition to Er. Of 
course this means that for some r values, one would be recovering the same Ej (i.e. those 
satisfying the condition \Ei\ < P, defined previously); however, since the EMM procedure 
is invariant under affine maps, this is just a redundancy. In this manner, bounds on r (or 
equivalently, Ej) and Eji would be generated. Our actual numerical results confirm this, 
and are consistent with the bounds quoted in the Tables. 



4/i8 + (34r 




3t 



(51) 
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We can enhance the above by including an additional constraint corresponding to 
Eq.(14). This is obtained as follows. ^Prom the Schrodinger equation, we know that 

f \^'(x)\^+i f {x^+ax)S(x) / .. , ^ ^ ■ ■ ^ f 4.U \ 

£/ = rkr\ ' (with respect to the imaginary part oi the energy) 

J ^(^) 

„ S{(i + Tr + a(i + T))S{i) 

= JW) • 

which reduces to (i.e. from Eq.(45)) 

//3 + 3Tfi2 + (3r^ + a)fj,i = 0. (53) 
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